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Atom- and site-resolved experiments with ultracold atoms in optical lattices provide a powerful 
platform for the simulation of strongly correlated materials. In this letter, we present a toolbox for 
the preparation, control and site-resolved detection of a tunnel-coupled hilayer degenerate quantum 
gas. Using a collisional blockade, we engineer occupation-dependent inter-plane transport which 
enables us to circumvent light-assisted pair loss during imaging and count n = 0 to n = 3 atoms 
per site. We obtain the first number- and site-resolved images of the Mott insulator “wedding cake” 
structure and observe the emergence of antiferromagnetic ordering across a magnetic quantum phase 
transition. We are further able to employ the bilayer system for spin-resolved readout of a mixture of 
two hyperfine states. This work opens the door to direct detection of entanglement and Kosterlitz- 
Thouless-type phase dynamics, as well as studies of coupled planar quantum materials. 


Introduetion 

Reduced and mixed dimensionality in solid state sys¬ 
tems is at the heart of exceptional material properties. 
Prominent examples include bilayer graphene mm , ex- 
citon condensation in bilayer systems [3] and uncon¬ 
ventional superconductors where superconductivity may 
originate from couplings in multilayer systems [4]. 

Experiments with ultracold atoms offer a clean and 
dissipation-free platform for the quantum simulation of 
such condensed matter systems [5]. Recently developed 
microscopy techniques with single-atom and single-site 
resolution provide direct access to local observables 
but so far have been constrained to two-dimensional sys¬ 
tems. Here, we present a scheme for high-fidelity fluores¬ 
cence imaging of a hilayer system with single-site resolu¬ 
tion. We realize full control over the resonantly coupled 
bilayer system and observe coherent dynamics between 
the two planes, confirming the suitability of our setup 
for investigations of strongly interacting bilayer materi¬ 
als. 

The bilayer system can be used to extend the ability of 
site-resolved optical lattice experiments to detect many- 
body ordering: Typically, the atomic hyperfine spin can¬ 
not be resolved during readout in quantum gas micro¬ 
scopes, leading to complications in the study of spin sys¬ 
tems H E]. More severely, only the parity of a lattice 
site occupation is accessible in fluorescence imaging due 
to light-assisted collisions and pairwise atom loss (“par¬ 
ity projection” ) IHlITlIIolIIi]. Several schemes for atom 
counting have been developed, providing global or aver¬ 
aged number statistics via spin-changing collisions [12] or 
an interaction blockade in double-wells m- We use this 
interaction blockade to engineer occupation-dependent 
transport between the two planes of the bilayer system. 
We circumvent parity projection and resolve lattice oc¬ 
cupation numbers n = 0 to n = 3 in one plane. Our 
technique allows for the first site-resolved, atom-number 
sensitive images of the Mott insulator “wedding cake” 
structure and of many-body ordering across a magnetic 
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FIG. 1. Preparation of a bilayer system. A degenerate gas 
of ®^Rb is loaded into two adjacent axial planes, near the fo¬ 
cus of a high-resolution imaging system. In combination with 
gravity, the axial lattice (spacing 1.5 /xm) and a 6x superlat- 
tice (spacing 9.2 jam.) result in a double-well geometry with 
tunable tunnel-coupling J and offset A. Both lattices along 
z are generated by reflecting beams off the flat surface of the 
hemispheric last lens of the imaging system. A square lattice 
projected through the objective onto the xy-plane realizes a 
bilayer Bose-Hubbard system. The optical molasses used for 
fluorescence imaging are also reflected at the surface, resulting 
in a standing wave pattern along z. 


quantum phase transition. Alternatively, we apply a 
magnetic field gradient to obtain spin-dependent trans¬ 
port between the planes and demonstrate spin-resolved 
readout in a mixture of two hyperfine spin states. 

The extension of quantum gas microscopy beyond the 
projection onto the atom number parity will give access 
to new observables. For one-and two-dimensional sys¬ 
tems, our techniques can yield complete number statistics 
and provide spectra of high-order density-density correla¬ 
tions and distribution functions. Such observables can be 
used to directly characterize quantum phases, for exam¬ 
ple the Tonks-Girardeau gas through its local and non¬ 
local pair correlations m or the dynamics of quantum 
phase transitions through counting statistics of excita¬ 
tions m- 
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Preparation of a resonant hilayer system 

The experimental setup, capable of single-atom re¬ 
solved in situ imaging of an ultracold cloud of ^^Rb, 
has been described in previous work m- Our experi¬ 
ments begin with a three-dimensional Bose-Einstein con¬ 
densate, which we compress in the direction of grav¬ 
ity (the ^-direction) in two-dimensional layers at the fo¬ 
cus of the imaging system. The confinement in the 2 ;- 
direction is provided by the “axial lattice,” with a spac¬ 
ing of d = 1.5 jam and a corresponding recoil energy of 
Er = 27 r X « 27r X 250 Hz. 

^ SrriRbd^ 

To bring two adjacent planes of the axial lattice into 
resonance, we employ a 6x superlattice with a spacing of 
9.2 jam. Both lattices are generated by reflecting beams 
off the flat surface of the last, hemispherical lens of the 
imaging system at large angles from normal (75° and 
87.6°). The relative phase between the two lattices at 
the position of the atoms can be tuned by changing the 
angle of incidence of the superlattice. In combination 
with the constant gradient of ^ = 27r x 3.2 g^xicd^lane 
gravity, we realize a resonant double-well system in the 
^-direction (Fig. [^. The residual offset A between the 
two axial planes of interest can be tuned by varying the 
depth of the superlattice, while the inter-plane tunnel¬ 
coupling J is controlled by the depth of the axial lattice. 
Other axial planes are sufficiently offset in energy that 
they remain entirely unpopulated. 

By controlling the initial position of the condensate, 
we can deterministically load one or two adjacent axial 
planes with arbitrary ratios of atom numbers. We use 
our imaging system to project a two-dimensional square 
lattice with spacing a = 680 nm onto the x^-plane, real¬ 
izing a bilayer Bose-Hubbard system. 

Imaging two planes with single-site resolution 

To image atoms in the bilayer system with single-site 
resolution, we use a sequential readout with two expo¬ 
sures. Separate high-resolution images of the two planes 
are obtained by shifting the focus of the imaging optics 
and tuning the fluorescence rate of atoms in different 
planes in real time. 

At the beginning of the imaging process, the atoms 
are localized and pinned in a deep optical lattice 60 GHz 
blue-detuned w.r.t the ^^Rb D1 line at 795 nm. They are 
simultaneously cooled by an optical molasses on the D2 
line (780 nm), and the scattered photons are collected 
to form images of the atom distribution. The molasses 
beams are reflected at an angle of 82° from normal off the 
flat surface of the hemisphere, forming a standing wave of 
period 2.8 jam along the ^-direction. This modulation of 
the molasses intensity corresponds to a modulation of flu¬ 
orescence rate for atoms in different planes. By changing 
the angle of incidence of the molasses with a galvanome- 
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FIG. 2. Site-resolved imaging of a bilayer quantum gas using 
two exposures. Image #1 (bottom row) comprises atoms in 
both planes, while image ff2 (top row) contains atoms only in 
plane I, directly revealing the corresponding atom positions. 
We obtain the atom distribution in plane II from image ffl 
after calculating and subtracting the background contributed 
by atoms in plane I. The procedure is illustrated here for de¬ 
coupled Mott insulators in both planes, with up to n = 2 
atoms per site. Framed panels denote the fitted atom distri¬ 
bution. 

ter, we are able to tune the ratio of fluorescence rates in 
planes I and H in real time. For the first exposure, we im¬ 
age both planes at a fluorescence ratio 1: 2 for 500 ms: 
Atoms in plane H with higher fluorescence rate are at 
the focus of the imaging system, while atoms in plane 
I contribute a weak out-of-focus background. Next, we 
remove atoms contained in plane H from the system. To 
this end, we increase the fluorescence ratio to 1: 3 and 
apply a higher molasses power for 300 ms. Due to the 
high scattering rate, atoms in plane H are now heated 
rather than cooled and ejected from the lattice, while 
atoms in plane I remained pinned. Simultaneously, we 
remove a 26 mm thick glass plate from the imaging path 
to shift the focus of the imaging system to plane I. At 
this point, plane H has been cleared of atoms and we take 
the second exposure to image atoms in plane I for 500 ms 
at an intermediate molasses power. 

The atom distributions in both planes are obtained 
after post-processing: We first determine the positions 
of atoms in the second image, containing only atoms 
in plane I. Subsequently we reconstruct the background 
contributed by these atoms to the first image by convolv¬ 
ing the extracted atom distribution with the measured 
point spread function of atoms one plane away from the 
focus. This simulated background is scaled to match the 
fluorescence rate of atoms in plane I during the first ex¬ 
posure and subtracted from the first image. A fit to the 
processed image then returns the atom positions in plane 
H. Figure illustrates the image analysis for decoupled 
Mott insulators in the two axial planes. 

The fidelity of the readout process is limited by our 
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ability to hold atoms in plane I while imaging and eject¬ 
ing atoms in plane 11. For optimized parameters, the life¬ 
time of atoms in plane I during the first exposure and 
ejection process is 27(2) s, resulting in a combined atom 
loss of 3.0% in plane I prior to imaging. The efficiency 
of the ejection of atoms from plane II is 99(1)%, lead¬ 
ing to an occasional unwanted background from atoms 
in plane II in the second image. In combination, these 
effects lead to an imaging fidelity of 95% in plane I and 
99% in plane II. The slight reduction in imaging fidelity 
for plane I primarily affects measurements in the Mott 
insulator phase, where quantum fluctuations of the lo¬ 
cal atom number are suppressed. Here, imaging errors 
effectively introduce an additional entropy of 0.2 ks per 
site, complicating access to physics at very low entropies, 
such as spin dynamics in the superexchange regime [8]. 
The current limitations on the imaging fidelity might be 
overcome by increasing the photon collection efficiency 
and choosing shorter exposure times, or by using addi¬ 
tional molasses beams on the D1 transition to decouple 
the cooling mechanism from the imaging process HZ]. 

Coherent dynamies between resonant axial planes 

We characterize the bilayer system by studying the 
double-well dynamics in the axial lattice direction. Our 
experiments begin with a single-layered Mott insulator 
in plane I in a deep two-dimensional optical lattice, ini¬ 
tially decoupled from plane II. Tunneling in the plane 
of the Mott insulator is negligible on time scales of our 
experiment and we concentrate on dynamics in the z- 
direction. 

We first set the depth of the superlattice to bring 
planes I and II near resonance. The axial lattice depth 
is then reduced from initially 250 Ej. to 8 in 2 ms, 
enabling inter-plane tunneling at a rate J 27r x 37 Hz. 
Figure a) shows offsets A at which particles may res¬ 
onantly tunnel from plane I to plane H. The first res¬ 
onance, at A — 27 r X 1.7 kHz corresponds to tunnel¬ 
ing from the ground band in plane I to the first excited 
band in plane H. The second resonance, at A = 0 corre¬ 
sponds to atoms tunneling within the ground band from 
plane I to plane H. For both processes, the on-site in¬ 
teraction shift of 1/ ^ 27 r X 300 Hz between singly and 
doubly occupied sites is well-resolved. Second-order tun¬ 
neling is expected to occur on a much slower timescale of 
^ ^ 27 r X 9 Hz and cannot be detected due to our parity 
projecting readout m- The horizontal scale is calibrated 
by measuring the offset A at various superlattice depths 
via photon-assisted tunneling [19]. 

We observe coherent Rabi oscillations between the two 
axial planes at the respective resonances for singly and 
doubly occupied sites, as shown in Fig. |^b). The ex¬ 
tracted resonant oscillation frequencies of the popula¬ 
tions are 118(7) Hz and 280(30) Hz respectively.The dy- 



FIG. 3. Inter-plane tunneling dynamics, a) Spectrum for 
inter-plane tunneling. Sites with single (n = 1, orange) 
and double (n = 2, blue) occupancy are initially prepared 
in plane I. After reducing the axial lattice depth to 7.9 Er 
(J ~ 27r X 37 Hz) in 2 ms and allowing atoms to tunnel for 
8 ms, we measure Podd, the probability to detect a single atom 
in plane I, as a function of offset A. Resonances correspond 
to tunneling into the ground band or first excited band of 
plane II, as indicated by sketches. The interaction shift be¬ 
tween n — 1 and n = 2 is clearly resolved (U ~ 27r x 300 Hz). 
Solid lines are Lorentzian fits to the data, b) Rabi oscilla¬ 
tions for n — 1 and n = 2 at their respective resonant offsets 
(A = 0 and A = 27 r x 300 Hz). Here, the axial lattice depth 
is ramped to 6.5 Er in 0.5 ms, giving J 27r x 55 Hz. Solid 
lines are damped sinusoidal fits from which the resonant Rabi 
frequencies are extracted. All error bars in this letter reflect 
1(7 statistical errors in the region-averaged mean podd- 

namics in singly occupied double-wells are in agreement 
with a numerical solution of the Schrodinger equation 
for single-particle wavefunctions and tunneling rates at 
an axial lattice depth of 6.5 while the rate of oscilla¬ 
tions on doubly occupied sites is enhanced by more than 
the factor \/2 expected from bosonic enhancement. We 
attribute this to the deformation of the potential when 
the tilt is applied using the superlattice, which can re¬ 
sult in reductions in both the effective double-well bar¬ 
rier height [20] and the effective double-well spacing. The 
observed damping of the Rabi oscillations is caused by 
inhomogeneities in the double-well potential across the 
two-dimensional plane. 


Beyond parity imaging 

Our technique of resonant population transfer to a sec¬ 
ond axial plane can be used to circumvent the limitations 
imposed by parity projection in optical lattice microscope 
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FIG. 4. Number-resolved observation of many-body ordering. 
After preparing a many-body state in plane I, occupation- 
sensitive transport of atoms to plane II allows the detection 
of occupancies n = 0 to n = 2. a) Processed single-shot image 
of the “wedding cake” structure of a two-shell Mott insulator, 
b) One-dimensional quantum phase transition from a para¬ 
magnetic phase to an antiferromagnetic phase. An increasing 
tilt is applied horizontally along three decoupled chains of 
length eight, tuning the system from unity filling (top) via 
the formation of doublon-hole pairs (middle) to a density- 
wave ordered state (bottom). 


experiments. We start by preparing a single-layered Mott 
insulator with singly and doubly occupied sites in plane I. 
With the axial tunnel coupling enabled (J ~ 27r x 48 Hz), 
we sweep the offset from A = 2.1 U to A = 0 in 75 ms, 
across the ground band tunneling resonance for doubly 
occupied sites at A « U. On doubly occupied sites, a 
single atom transitions at an offset corresponding to the 
on-site interaction U. The transfer of a second atom 
is suppressed by a collisional interaction blockade, leav¬ 
ing one atom in plane I and one atom in plane II [13]. 
Atoms on singly occupied sites distribute over planes I 
and II with roughly equal probabilities. At the end of 
the sweep we image both planes and obtain the distri¬ 
bution of holes, single atoms and doublons in the initial 
Mott insulator by adding the atom distributions from 
both planes. The reconstructed “wedding cake” struc¬ 
ture of a Mott insulator is shown in Fig. |^a), combining 
single-site resolution igiT] and atom-number sensitive de¬ 
tection [m . Taking into account the defect probabil¬ 
ity of the initial Mott insulator and the bilayer imaging 
fidelity, we achieve 96(1)% efficiency of separating the 
doublons into two planes. 

We employ our imaging technique to detect many-body 
ordering across a magnetic quantum phase transition de¬ 
scribed in previous work |20l [23] . Our experiments start 
with a n = 1 Mott insulator in plane I, decoupled into 
one-dimensional chains along the x-direction. Using a 
magnetic field gradient along the chains , we drive a quan¬ 
tum phase transition from a paramagnetic state (unity 
filling) to an antiferromagnetic (density-wave ordered) 
state. We image the atom distribution at various points 
along the transition, carrying out the beyond-parity read- 
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FIG. 5. Resolving up to three atoms per site in a binary read¬ 
out scheme, a) Occupations n = 1 to n = 3 in plane I are 
mapped to different distributions in plane I and II by inter¬ 
plane transfer and subsequent parity projection, b) Averaged 
Podd in plane I (blue) and plane II (green) after preparing a 
three-shell Mott insulator in plane I and mapping the occupa¬ 
tion onto the two planes. Mott-insulating regions from n = 3 
at the trap center to n = 1 near the trap perimeter are re¬ 
solved. Gurves are fits with (concatenated) error functions. 


out scheme as above. The formation of doublon-hole 
pairs and antiferromagnetic ordering is visible in single¬ 
shot reconstructions of the atom distribution in Fig.j^b). 
In contrast to the previous global detection of antifer¬ 
romagnetic order [23], the ability to resolve individual 
doublon-hole pairs enables direct measurement of the 
Neel order parameter, and detailed studies of phenom¬ 
ena such as frustration and the dynamics of defect or 
domain formation in the underlying model. 

A further generalization of our readout scheme allows 
the unambiguous detection of atom numbers n = 0 to 
n = 3. Using each side of the double-well as a “bit” 
that is either bright (odd occupancy) or dark (even oc¬ 
cupancy) after parity projection, four different number 
states can be encoded. Figure]^ a) illustrates the map¬ 
ping for this “binary readout”. After preparing a three- 
shell Mott insulator in plane I, we ramp the offset A 
through the resonances for n = 2 and n = 3 atoms per 
site, stopping before reaching the n = 1 resonance at 
A = 0. Figure l^b) shows Podd after the ramp vs. radial 
distance for both planes. All plateaus of constant atom 
number from n = 3 at the center of the cloud to n = 1 
on the outside edge can be identified. 

The fidelity of the atom-number sensitive readout is 
























































5 


currently limited by the small energy scales for dynam¬ 
ics in the ^-direction. The relatively small interaction 
{U ^ 27r X 300 Hz) and the large spacing of the axial lat¬ 
tice (1.5 yum) lead to slow dynamics and sensitivity to 
lattice inhomogeneity. By using a Feshbach resonance 
and a smaller axial lattice spacing, the robustness of the 
mapping process onto axial planes could be further im¬ 
proved. 


Spin-resolved readout 




plane © (j) 


b 



The presence of a spin degree of freedom greatly ex¬ 
tends the capability of experiments with ultracold atoms 
to simulate condensed matter Hamiltonians. Our bi¬ 
layer system may be used to observe spin ordering in 
such systems, providing spin-sensitive readout in a two- 
dimensional mixture of two spin states. The scheme is 
illustrated in Fig.[^ A mixture of two appropriately cho¬ 
sen hyperfine states initially resides in plane I. To map 
out the distribution of both spin states in plane I, we 
enable transport between the two planes, after motion 
within the planes has been frozen out by a deep lattice. A 
magnetic field gradient applied in the ^-direction causes 
atoms in one hyperfine state to transfer to plane H, while 
atoms in the second hyperfine state experience a force in 
the opposite direction and remain in plane L The hyper¬ 
fine spin degree of freedom is thus mapped to the two 
planes of the axial lattice, and both spin states can be 
imaged simultaneously. 

We demonstrate spin-resolved readout for the two hy¬ 
perfine states \F^mF) = If,—1) (labeled |t)) and |2,—2) 
(ID), for which and Qf'^f = —1, respec¬ 

tively. Initially, we prepare a n = 1 Mott insulator of 
atoms in |t) in plane I, and transfer atoms to the |D 
state using a resonant microwave pulse in a bias field of 
1.5 G. After adiabatically transferring the atoms from 
the axial lattice into a single well of the superlattice, we 
ramp up a magnetic field gradient of 30 G/cm in the 2 ;- 
direction in 70 ms. When the axial lattice is then ramped 
back on, the magnetic field gradient causes atoms in state 
1^) to transfer to plane H, while atoms in state |t) re¬ 
main in plane I. Figure shows the population of both 
planes after mapping versus microwave pulse duration. 
The sinusoidal variation in anti-phase demonstrates the 
mapping of spin to plane degree of freedom. We estimate 
a fidelity of 97(1)% for the correct sorting of hyperfine 
spin into different axial lattice planes from the offset and 
amplitude of the fit. 

Unlike previous experiments, in which only one of two 
spin states could be imaged in situ, our technique gives 
access to the full spin distribution in an interacting many- 
body system. This scheme will enable further studies of 
two-component systems, such as impurity dynamics [8] 
and collective excitations m- 


FIG. 6. Spin-resolved readout, a) Procedure for mapping 
the hyperfine spin onto different axial planes: After reducing 
the tunnel barrier in a double-well with an arbitrary spin in 
plane I to zero, a B-field gradient separates the |t) and jj.) 
components into different axial planes, b) A resonant mi¬ 
crowave pulse is applied to a n = 1 Mott insulator in state 
It) in plane I. The sinusoidal variation in hyperfine spin is 
mapped onto occupation in plane I (blue) and plane II (green). 
Curves are fits with sine functions. 

Summary 

We have demonstrated single-atom and single-site re¬ 
solved imaging of a bilayer quantum degenerate gas. We 
engineer occupation-sensitive transport between resonant 
layers with an optical superlattice and the use of a col- 
lisional interaction blockade. This transport and read¬ 
out scheme circumvents the problem of parity projec¬ 
tion in quantum gas microscopes, allowing the unam¬ 
biguous identification of atom numbers up to n = 3. 
We have obtained the first number- and site-resolved 
images of the Mott-insulator “wedding-cake” structure 
and directly imaged many-body ordering across a quan¬ 
tum phase transition. Making use of a magnetic field 
gradient in the ^-direction, we have demonstrated spin- 
dependent transport and spin- and site-resolved readout 
of a two-species mixture. This approach will facilitate the 
observation of antiferromagnetic ordering in the Fermi- 
Hubbard model and spin-dependent phenomena such as 
spin-charge separation m- 

Direct extensions to our imaging scheme will enable 
site-resolved detection in more than two planes. In par¬ 
ticular, the use of additional molasses beams on the 
D1 transition for cooling m while collecting fluorescence 
photons only from molasses on the D2 transition will lead 
to even more spatially selective readout. 

Our techniques for the preparation and readout of res¬ 
onant bilayer systems open numerous possibilities for fur¬ 
ther studies of low-dimensional phenomena: Interfering 
two planar superfluids should enable in situ observation 
of phase evolution in two dimensions and the dynam¬ 
ics of the Kosterlitz-Thouless phase transition [25] , while 
many-body entanglement can be measured in a system of 
two copies of a planar system m- Finally, our technique 
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could be used to reduce the entropy in two-dimensional 
Mott insulators by “filling” empty sites (defects) with 
atoms by merging with a reservoir plane m 

We thank Andrew Daley for helpful discussions. This 
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Note added: Spin-resolved detection of individual 
atoms in a single one-dimensional Hubbard chain was 
recently reported in Ref. [28] 
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